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PROCEDE DE RECONSTRUCTION D'UNE IRflAGE TRIDIMENSIONNELLE D'UN OBJET. EN PARTICULIER 
UNE IMAGE TRIDIMENSIONNELLE ANGIOGRAPHIQUE. 

Le proc6d6 comprend un 6talonnage de rappareil 
dans lequel on 6labore un volume virtue! entourant I'objet et 
d6compos6 en voxels Vi, une acquisition dudit ensemble 
d'images projet^es bidimensionnelles num6ris6es, et une 
reconstruction 33 de I'lmage tridimensionnelle ^ partir des 
images bidimensionnelles projet^es acquises et d'un algo- 
rithme iteratif de reconstruction alg^brique d'image. On ef- 
fectue une premiere iteration de I'algorithme avec une 
resolution d'image initiale pr6d6termin6e de facon k obtenir 
k Tissue de cette premiere iteration des premieres valeurs 
de densites pour les voxels dudit volume Wl , on subdivise 
44 une partie au moins W2 desdits voxels du volume virtuel 
en plusieurs ensembles VV3, W5, W6, W7 correspon- 
dant respectivement k des resolutions d'images diff6rentes 
multiples ou sous-multiples de la resolution initiale, et au 
sein de chaque iteration suivante de Talgorithme, on appli- 
que successtvement ledit atgorithme k partir de chacun des- 
dits ensembles de voxels. 
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Precede de reconstruction d*une image tridimensionnelle d*un 
objet, en particulser une image tridimensionnelle angiographique* 

L'invention conceme la reconstruction d'une image 
tridimensionnelle d'un objet k partir d'un ensemble d'images projet^es 
bidimensionnelles dudit objet obtenues pour diff6rentes positions d'un 
appareil de prises de vues autour de Tobjet. 

Elle trpuve une application particuli^rement int^ressante dans le 
domaine medical dans lequel on procfede k la reconstruction des structures 
intemes de patient sous examen, en particulier la reconstruction d'images 
angiographiques, c*est-i-dire Fobtention d'images d'arbres vasculaires 
opacifies par injection d'un produit de contraste. 

L'inv^entio^^ ^ ?^'9^y!^r_n6anr^^ dans 

d'autres domaines, notamment dans le contrdle non destructif industriel 

dans lequel des examens du mSme type que lesexamens^ 

pratiques. 

Dans le domaine medical, les images projet^es 
bidimensionnelles de Tobjet, par exemple la tfite d'un patient, sont 
gen6ralement obtenues par la rotation d'un appareil de prises de vues k 
rayons X toumant autour de I'objet. 

II existe essentiellement deux types d'algorithmes de 
reconstruction en imagerie par rayons X. 

Un premier type pr6voit un calcul de r^troprojection et de 
filtrage ou encore une reconstruction par transform^e de Fourier k 
plusieurs dimensions. 

Un deuxi^me type, celui qui est concern^ par I'invention, 
conceme les m^thodes it6ratives de reconstruction dites encore 
algdbriques. Le principe d'un tel algorithme alg^brique est parfaitement 
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connu de Thomme du metier et a fait d6jk Tobjet de nombreuses 
publications. On citera notamment Tarticle de GORDON, bender, et 
HERMAN intitul6 "Algebraic reconstruction technic for three dimensional 
Electron nucrocopy and X ray photography". Journal THEO. BIOL. 29, 
pages 471 ^ 781 (1970), ou bien Touvrage de Anil K. JAIN intitule 
"Fondamentals of digital image processing" Prentice Hall information 
and system sciences series, Thomas Kailath Series Edition, ou encore les 
demandes de brevets fran^ais n° 89 03606 ou n^ 89 16906. 

En bref, aprfes un etalonnage de Tappareil utilis6 pour 
determiner notamment les param&tres de la projection dans les plans de 
projection des images acquises, d'un volume observe d6compos6 en 
616ments volumiques 61^mentaires ou voxels (ces paramfetres 
d'6talonnage formant des matrices de projection), Talgorithme de 
reconstruction alg6brique d'image est utilise pour reconstruire le volume 
tridimensionnel k partir de ces images projet6es bidimensionnelles. Le 
principe de base de cet algorithme est d*initialiser les voxels du volume ^ 
une valeur initiale pr6d6termin6e, par exemple une valeur nulle, et d'it^rer 
un certain nombre de fois les operations suivantes : projection des voxels 
dans le plan de chaque image acquise de fa9on k obtenir une image 
virtuelle, determination de la difference entre le volume projet6 (image 
virtuelle) et Timage acquise correspondante puis retroprojection de cette 
difference dans le^ volume. Aprfes un certain nombre^ d*it6rations,. on 
obtient pour chaque voxel une valeur estimde representative de la density 
de produit de contraste inject6 dans les vaisseaux radiographics ce qui 
permet de visualiser en trois dimensions la cartographic de ces vaisseaux 
radiographics. 

Les images acquises ont gCnCralement une resolution Cgale k 
512, c'est-2l-dire qu'elles comportent 512 lignes et 512 colonnes de pixels. 
Si Ton appliquait Talgorithme de reconstruction d*images sur la totalite 
des images acquises, on obtiendrait un traitement s'effectuant sur environ 
128 nnillions de voxels, ce qui est un nombre beaucoup trop eiev6 et de 
toutes fa^ons pas tr&s utile car les structures vasculaires que Ton souhaite 
genCralement visualiser occupent typiquement environ 2% du volume 
virtuel, 

Aussi, il a ete propose de diminuer la valeur de resolution des 
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images acquises de fafon, en effectuant des moyennes de quatre pixels, k 
se ramener k une resolution de valeur 256, ce qui conduit k une reduction 
du volume virtuel. Apr&s une premifere iteration de Talgorithme effectu6e 
sur chaque image, on obtient d6j^ une premiere representation grossifere 
de la localisation des objets d*int6ret, ce qui permet de s61ectionner, pour 
les iterations suivantes, un sous-ensemble de p voxels et d'61iminer les 
autres. Chacun de ces p voxels restants sont alors subdivis6s en 8, de fagon 
k ce que la projection d'un tel voxel subdivise corresponde k 1/2 pixel de 
resolution 256, c'est-k-dire k un pixel de resolution 512. On effectue alors 
les iterations suivantes de I'algorithme, en pratique deux iterations, sur 
Tensemble de ces voxels subdi vises. 

Or, meme aprfes reiimination d'un certain nombre de voxels, 
Tensemble des voxels subdivis6s se monte k environ 32 millions, ce qui 
reste un nombre trfes important, et se traduit par un codt non negligeable en 
temps de calcul. Une solution pourrait consister alors k eiiminer plus de 
voxels pour diminuer encore le nombre de voxels subdivises restants. 
Mais, si Ton va dans cette voie, on obtient alors des artefacts sur les 
images visualisees, en raison de la sous-representation du volume des 
donnees prises en compte par Talgorithme de reconstruction d'image. 

L'inventiqn^^ k apport er une solution ^ ce probldme. 

Un but de Tinvention est de reduire le temps de calcul du 
processeur impiement^^ d'image tout en 

ne modifiant pas la qualite des images obtenues, c*est-^-dire en 
n*introduisant pas notanunent d'artefacts dans les images. 

L'invention propose done un procede de reconstruction d une 
image tridimensionnelle d'un objet k partir d'un ensemble d'images 
projetees bidimensionnelles numerisees dudit objet obtenues par 
differentes positions d*un appareil de prises de vues autour de Tobjet. Ce 
procede comprend un etalonnage de I'appareil dans lequel on eiabore un 
volume virtuel entourant I'objet et decompose en voxels, une acquisition 
dudit ensemble damages projetees bidimensionnelles numerisees, et une 
reconstruction de I'image tridimensionnelle k partir des images 
bidimensionnelles projetees acquises et d'un algorithme iteratif de 
reconstruction algebrique d'image. 

Selon une caracteristique generale de Tinvention, on effectue 
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une premifere iteration de ralgorithme avec une resolution d*image initiale 
pr6d6termin6e de fa^on k obtenir h Tissue de cette premifere iteration des 
premieres valeurs de density pour les voxels dudit volume. On subdivise 
une partie au moins desdits voxels du volume virtuel en plusieurs 
ensembles correspondant respectivement h des resolutions d'image 
diffSrentes, multiples ou sous-multiples de la resolution initiale. Et, au 
sein de chaque iteration suivante de Talgorithme, on applique 
successivement ledit algorithme h partir de chacun desdits ensembles de 
voxels. 

En d'autres termes, Tinvention pr^voit d'appliquer Talgorithme 
iteratif de reconstruction alg^brique d'images sur un volume multi- 
resolutions. Dans ce volume, on seiectionne des voxels qui sont a priori 
representatifs des objets d'interet h visualiser et on les divise de fa9on h 
augmenter la resolution. Les autres voxels, qui sont moins interessants, 
car ne concemant pas directement les objets h visualiser, sont soit laisses 
identiques, soit regroupes au moins une fois de fagon k diminuer la valeur 
de la resolution, mais sont toutefois utilises pour les calculs de 
reconstitution d'image, ce qui permet d'obtenir en final des images de trhs 
bonne qualite avec un temps de calcul reduit. 

Selon un 4?Jl^y®??i^^%_^^ cours de 

chaque iteration suivante de Falgorithme, c'est-^-dire au cours de chaque 
iteration k partir de la deuxifeme, on determine pour chaque image projetee 
acquise, une image virtuelle en sommant successivement des images 
virtuelles eiementaires respectivement obtenues par les projections dans 
le plan d'acquisition de ladite image projetee acquise des ensembles 
correspondants de voxels. 

Plus precisement, et notanunent dans le but de reduire encore le 
temps de calcul, on eiabore de preference les images virtuelles 
eiementaires dans Tordre croissant des valeurs de resolution en 
conmien9ant par Timage virtuelle eiementaire correspondant k la valeur 
de resolution la plus faible (par exemple 64), et, aprfes avoir determine une 
image virtuelle eiementaire, on effectue une modiflcation d'echelle de 
fa9on k obtenir une image virtuelle dite "agrandie" dont la resolution 
correspond k celle de la prochaine image virtuelle eiementaire k 
determiner, et on somme ladite image virtuelle agrandie avec ladite 
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prochaine image virtuelle 616mentaire, 

D'une fa^on g6n6rale, selon un mode de mise en oeuvre de 
rinvention, on acquiert des images projet6es bidimensionnelles ayant une 
resolution de base pred6temun6e (par exemple r = 512). On choisit une 
5 resolution initiale 6gale k un sous-multiple de la resolution de base (par 
exemple r/2 = 256), et on choisit parmi lesdites resolutions differentes 
d'image, ladite resolution de base (r), la resolution initiale (r/2) et au 
moins une premiere resolution suppiementaire (par exemple r/4 = 128) 
sous-multiple de la resolution initiale (r/2). 

10 Selon un mode de mise en oeuvre de Tinvention, on eiabore un 

premier seuil de densite en fonction d'un critfere de selection 
predetermine. On subdivise chaque voxel ayant une densite superieure ou 
egale au premier seuil, en un premier nombre de voxels subdivises, le 
premier nombre (typiquement 8) etant defini compte tenu du rapport entre 

15 la resolution de base et la resolution initiale, tous les voxels subdivises 
formant un premier ensemble de voxels correspondant h ladite resolution 
de base (r = 5 1 2). On regroupe certains au moins des voxels dont la densite 
est inferieure audit premier seuil et qui satisfont k un critfere de 
regroupement predetermine, de fa9on k former des voxels regroupes 

20 formant un deuxieme ensemble de voxels correspondant k ladite premiere 

resolution suppiementaire (r/4 = 128), le nombre (typiquement 8) de 
voxels regroupesf dans chaque groupe tenant compte du rapport entre la 
resolution initiale (r = 256) et ladite premifere resolution suppiementaire 
(128). Et, les voxels dont la densite est inferieure audit premier seuil et qui 

25 ne satisfont pas au entire de regroupement predetermine forment un 
troisieme ensemble de voxels correspondant k ladite resolution initiale 
(256). 

A titre d'exemple, un voxel satisfait audit critfere de 
regroupement si chacune des coordonnees du centre du voxel est un 
30 multiple de 2 et si la densite de chaque voxel voisin dudit voxel est 
inferieure audit premier seuil. 

On peut egalement, de preference, regrouper une fois de plus les 
voxels d6jk regroupes, et qui satisfont audit critfere de regroupement, de 
fa5on k former un quatrifeme ensemble de voxels correspondant k une 
35 deuxifeme resolution suppiementaire (par exemple 64) sous-multiple de la 
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premiere resolution suppl^mentaire (par exemple 128). 

D*une fa^on g6n6rale, la partie des voxels du volume virtuel, que 
Ton subdivise, peut etre obtenue en 61iminant les voxels situ6s dans une 
couche d'^paisseur pr6d€tennin6e de la surface du volume. II a en effet 6t6 
observe que Talgorithme de reconstruction d'image tendait k cr6er des 
fortes densit^s sur le pourtour du volume et qui en fait ne correspondent 
pas k des objets d'intirgt k visualiser. 

On peut egalement eiiminer du volume virtuel, les voxels isol6s 
dont la density est sup6rieure k un deuxi&me seuil predetermine. En effet, 
il a ete observe qu'un voxel isoie de forte densite ne correspondeiit pas k un 
objet d'interSt. 

Par ailleurs, aprfes la premiere iteration de Talgorithme, on 
affecte de preference k chaque voxel dont la valeur de densite est 
superieure ou egale k un troisifeme seuil predetermine, une valeur de 
densite egale k la valeur maximale de densite prise parmi les valeurs de 
densite dudit voxel et des voxels voisins. On effectue en d'autres termes 
ici, de preference, une dilatation des voxels de forte densite, de fa9on k 
ameiiorer encore la qualite de la visualisation des objets d'interet. 

D'autres avantages et caracteristiques de I'invention 
apparaitront ^ rexamen de la des^^^^ detailiee d'un mode de noise en 
oeuvre nuUement limitatif, et des dessins annexes, sur lesquels : 

- la figure 1 illustre schematiquement un jeu d'images projetees 

bidimensionnelles autour d'un objet, 

- la figure 2 illustre plus en detail Tacquisition de Tune de ces 
images projetees bidimensionnelles, 

- la figure 3 est un organigranmie d'un mode de mise en oeuvre du 
procede selon Tinvention, et 

- la figure 4 illustre plus en detail une partie de Torganigramme 
de la figure 3. 

Bien que Tinvention n'y soit pas limitee, on va maintenant 
decrire I'application du procede selon I'invention k la reconstruction d'une 
image angiographique tridimensionnelle d'un patient, en particulier sa 
tete. 

Si Ton se reffere plus particuli^rement aux figures 1 et 2, on voit 
que le systfeme d'imagerie utilisable pour mettre en oeuvre I'invention 
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permet d'obtenir un ensemble dlmages acquises bidimensionnelles lAl- 
lAn obtenues en Tepfece par la rotation autour de la tete 1 d'un patient d'une 
source 2 k rayons X. En fait, comme il est classique en angiographic, 
chaque image acquise LAi est une image soustraite qui est par exemple 
obtenue par une technique classique de soustraction logarithmique de 
deux radiographics prises sous la m§me incidence avant et aprfes une 
injection d'un produit de contraste dans Tarbre vasculaire dont on souhaite 
reconstruire Timage tridimensionnelle. 

Chaque image acquise lAi est obtenue k partir d'un detecteur de 
rayonnement bidimensionnel, par exemple du type ampliftcateur de 
luminance utilise en radiologic, dispos6 en vis-k-vis du tube h rayons X 
dans un plan dit plan de projection PPi. Les diffdrents plans de projection 
sont obtenus par les diff6rentes positions angulaires du detecteur en 
rotation autour de la tSte du patient. La normale XOi au plan de projection 
PPi d^finit I'axe optique de I'image acquise lAi. Le detecteur est relie h des 
moyens de traitement 3 comportant notamment des moyens 
d'echantillonnage relics k un microprocesseur incorporant de fa9on 
logicielle dans sa mdmoire de programme associ^ I'algorithme de 
reconstruction algebrique d'images utilise dans Tinvention et d'une fagon 
g6n6rale tous^les moyens fonctiojinels permettant la mise en oeuvre du 
proc6d6 scion I'invention. 

Dans le cas d'un systfeme d'imagerie^^ rayons X compost d'une 
source de rayons X et d'un detecteur bidimensionnel, I'operation 
g6om6trique intervenant dans la production de I'image acquise est une 
projection conique d'un objet analyst, d6ploy6 dans un espace k trois 
dimensions, sur un espace k deux dimensions qui est celui du plan de 
projection correspondant au plan de detection. Les param&tres 
g6om6triques d6crivant compl^tement les divcrses projections coniques 
doivent Stre connus. Or, il est souvent impossible et trop impr6cis, 
d'acc^der k ces paramfetres de fa9on directe, c'est-^-dire, par exemple, en 
mesurant directement sur le systfeme d'acquisition la distance entre la 
source de rayons X et le detecteur. 

On appelle "6talonnage" d'un systfeme d'imagerie, I'operation 
qui aboutit k la connaissance indirccte precise des param&tres 
g^om^triqucs qui intcrviennent dans la production d'une image. Le 
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principe, classique et connu, est bas6 sur rutilisation d'un fantdme 
geom6trique connu dans Tespace tridimensionnel, et dont on acquiert la 
projection bidimensionnelle. Plus pr^cis6ment, I'etalonnage comporte les 
6tapes suivantes : 

- on dispose d*un objet connu, le fantSme d*6talonnage, 
pr^sentant un certain nombre de points caract6ristique$ dont la position 
dans Tespace est connue par des coordonndes mesur6es par rapport k un 
repfere propre k cet objet; 

- on acquiert Timage de ce fantdme dans les conditions 
g6om6triques d'un point de vue (ou incidence) que Ton veut 6talonner; 

- on reconnait les projections des points caract^ristiques dans 
rimage. Pour cela on associe chaque point caract^ristique de Tobjet h sa 
trace dans Timage acquise projet^e; 

- on inverse, au sens math6matique, le systfeme d*6quation 
d^crivant la projection; 

- et on obtient finalement Tensemble des param^tres de la 
projection pour le point de vue donne. 

Une forme de fantdme d'6talonnage geom^trique souvent 
utilis^e est celle d'un cube, aux huit coins duquel sont dispos6es des billes 
m^talUques opaflue^ jaux rayons >C JL'^taloM operation 
connue de Thonune de metier, elle ne sera pas d6crite plus en details, 
_d*autant que plusieurs publications, ont d6j^ d^crit le principe d'une 
calibration g6om6trique manuelle. On peut notamment citer les articles 
suivant : 

(1) D.L. Parker, J.Wu, D.L Pope, R. Van Bree, G.R. Caputp and 
H.W. Marshall, "Three-dimensional reconstruction and flow 
measurements of coronary arteries using multi view digital angiography", 
in New Developmnts in Quantitative Coronary Arterography , J.C. Reiber 
and P.W. Serruys Eds., pp. 225-247, Kluwer Academic Publishers, 1988; 

(2) D.J. Hawks, A.C.F. Colchester and C.R. Mol; "The accurate 3-D 
reconstruction of the geometric configuration of the vascular trees from 
X-ray recordings," in Physics and Engineering of Medical Imaging, 
R. Guzzardi Ed., Nijhoff, 1987. 

(3) M. Garreau, J-L Coatrieux, R. Collorec and C. Chardenon, 
"A Knowledge-based approach for 3-d reconstruction and labeling of 
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vascular networks from biplane angiographic projections", EEEE Medical 
Imaging, vol. 10, n^2, pp. 122-131, June 1991. 

II est 6galement possible d'utiliser un proc6d6 connu 
d'^talonnage g6om6trique automatique d*un syst^me d'imagerie par 
rayons X, tel que celui d6crit dans la demande de brevet fran9ais 
n° 93 00804. En bref, on utilise pour un tel €talonnage automatique un 
fantdme dans lequel les billes sont distributes, de proche en proche, en 
une succession telle que des altitudes de billes, mesurtes le long de Taxe 
de rotation du systfeme d'imagerie, et surtout un axe du fantome, soieht 
monotones croissantes (ou d6croissantes) avec un num6ro d'ordre des 
billes dans la succession. 

L'ttalonnage du systfeme d'imagerie permet notamment de 
determiner I'axe moyen estimt Ax de rotation de I'appareil de prises de 
vues autour de I'objet ainsi que la position de la source 2 et les 
caracttristiques gdom^triques des axes optiques des diffSrentes images 
acquises. L'6talonnage permet 6galement de dtfinir un volume virtuel VV 
(intersection des difftrents c6nes de projection) entourant Tobjet 1 et 
d6compos6 en elements volumiques 616mentaires Vi ou "voxels". Ce 
volume VV, et done chaque voxel Vi, est spatialement rep6re dans un 
ref6rentiel, dtnomnie ci-aprfes r6f6rentiel de calibration, dont J'un des 
axes, en I'espfece Taxe Z, est confondu avec Taxe de rotation estim6 Ax. II 
convient de noter ici que les plans de projection PPi dans lesquels se 
projettent les images acquises lAi ne sont g6n6ralement pas parall&les h 
Taxe Z. 

L'etalonnage permet 6galement de d6finir pour chaque image 
acquise lAi une matrice de projection Pi permettant de dttemruner, pour 
chaque voxel Vi les coordonn6es de sa projection (pixel) dans Timage 
acquise lAi correspondante. 

On va maintenant dtcrire en se r6f6rant plus particuli&rement 
aux figures 3 et 4, un mode de mise en oeuvre du proc6d6 selon I'invention. 

L'algorithme de reconstruction algtbrique d'image est appliqu6 
directement sur les images acquises lAi (obtenues aprbs ttalonnage 30 et 
acquisition 31) de fa9on k d61ivrer I'image tridimensionnelle reconstitute 
IF, typiquement apr&s trois iterations. 

On rappelle ici que le principe de base de l'algorithme de 
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reconstruction d'image selon Tinvention consiste k initialiser les voxels 
du volume k une valeur initiale pr6d6tennin6e, par exemple la valeur 
nulle, et d'effectuer ensuite un certain nombre dlterations, typiquement 
trois. Chaque iteration comporte, pour chaque image acquise, la 
projection des voxels dans le plan de chaque image acquise de fa^on k 
obtenir une image virtuelle, la determination de la difference entre le 
volume projet6 (image virtuelle) et I'image acquise correspondante, puis 
la retroprojection de cette difference dans le volume. 

G6n6ralement aprSs trois iterations, on obtient pour chaque 
voxel une valeur estim^e representative de la densite de produit de 
contraste injecte dans les vaisseaux radiographies, ce qui permet de 
visualiser en trois dimensions la cartographic de ces vaisseaux 
radiographies (image IF). 

Comme on le voit sur la figure 4, alors que la resolution de base 
des images acquises est typiquement egale ^512(512x512 pixels), on 
ramfene cette resolution de base k une resolution initiale moitie (256), ce 
qui revient k regrouper les voxels du volume par huit et k effectuer sur 
chaque image acquise un regroupement de pixels par groupes de quatre. 

On effectue alors une premifere iteration (etape 40) et Ton obtient 
pour chacun des voxels du volume VVl, des premieres valeurs de densite. 

Alors que Ton pourrait en theorie appliquer la deuxieme 
iteration directement sur ce premier volume VVl, il est particulierement 
avantageux, de fa^on k diminuer le temps de calcul, d'eiiminer dej^ une 
partie des voxels de ce premier volume VVl. 

Plus precisement, on procSde k une elimination des voxels du 
pourtour du volume VVl (etape 41). En pratique, cela revient k "peler" le 
volume, c'est-^-dire k retirer une couche superficielle de la surface du 
volume ayant par exemple une "epaisseur" de 3 voxels. 

On procfede ensuite k une elimination des voxels brillants isoies 
(etape 42). En d'autres termes, on eiimine tous les voxels ayant une valeur 
de densite superieure k un seuil predetermine et dont les voxels voisins ne 
presentent pas des valeurs de densite superieures k ce seuil. 

On effectue ensuite (etape 43) une dilatation des voxels 
brillants, c'est-^-dire une dilatation des voxels dont la valeur de densite 
est superieure k un seuil predetermine t3. En pratique, on affecte k chaque 
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voxel dont la valeur de density est superieure au seuil 13, la valeur 
maximale de density prise parmi la valeur de density dudit voxel et les 
valeurs de density des voxels voisins de ce voxel. 

On obtient, aprfes cette phase de pretraitement, un volume restant 
5 de voxels VV2 sur lequel on va effectuer deux iterations successives de 
Talgorithme de reconstruction d'image. 

Avant d'effectuer la deuxifeme iteration, on s61ectionne parmi 
les voxels du volume V V2, ceux dont la valeur de density est sup6rieure k 
un premier seuil predetermine tp On subdivise alors ensuite ces voxels en 

10 8 (etape 44) de fajon k obtenir des voxels subdivisfis formant un premier 
sous-ensemble VV3 et correspondant k une resolution d'image double de 
la resolution initiale, c'est-k-dire k la resolution de base egale k 512. 
Parmi les voxels non subdivis6s du volume VV2, on seiectionne ceux qui 
peuvent etre regroupes compte tenu d*un entire de regroupement (etape 

15 45). Get ensemble de voxels regroupes par groupes de huit, va etre destine 
k former un deuxifeme sous-ensemble de voxels correspondant k une 
resolution moiti6 de la resolution initiale, c'est-i-dire k une resolution de 
128, 

Le critfere de regroupement est le suivant : 

20 _ . ^ peut etre regroupe avec les sept voxels voisins si les 

coordonnees de son centre sont chacune un multiple de deux et si la valeur 
de la densite de chacun des sept voisins du voxel est egalement inferieure 
au seuil tp Si Tune ou Tautre de ces deux conditions n'est pas realisee, 
alors meme si un voxel a une densite inferieure au seuil tp il ne peut pas 

25 Stre regroupe et restera done tel quel. A Tissue de ce premier 
regroupement, on obtient done un sous-ensemble VV4 de voxels 
regroupes correspondant k une resolution de 128, ainsi qu'un sous- 
ensemble VV7 de voxels non regroupes et ayant une densite inferieure au 
premier seuil tj et correspondant done k une resolution de 256. 

30 On va ensuite appliquer avantageusement le mSme critfere de 

regroupement au sous-ensemble VV4, de fagon k effectuer un deuxifeme 
regroupement de voxels d6jk regroupes, de fa^on k obtenir un sous- 
ensemble VV6 de voxels issus du sous-ensemble VV4 et regroupes encore 
par groupes de huit, et correspondant done k une resolution de 64. 

35 Puisque le mSme critfere de regroupement est applique, un voxel 
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du sous-ensemble VV4 pouira Stre regroup^ avec ses sept voisins si ses 
coordonnees sont un multiple de deux et si chacun de ses sept voisins a 
egalement une densit6 inf6rieure au premier seuil tj. 

Les voxels du sous-ensemble VV4 qui n'ont pas pu §tre 
regroup^s une deuxifeme fois forment alors le sous-ensemble VV5 
correspondant done k la resolution 128. 

On procfede ensuite (6tape 47) h T^laboration, pour chaque image 
acquise, d'une image virtuelle correspondant h la projection des voxels 
dans le plan d'acquisition de ladite image acquise. 

Ceci ^tant, etant donn6 que le volume sur lequel on effectue cette 
projection est un volume multi-r6solutions, on va en fait effectuer des 
projections successives k partir des sous-ensembles de voxels successifs 
de fa9on k ddfinir des images virtuelles 616mentaires que Ton va agrandir 
puis sommer aux projections 616mentaires correspondant aux sous- 
ensembles suivants de voxels. 

Plus pr6cis6ment, k partir du sous-ensemble V V6 correspondant 
k la resolution 64, on eiabore par projection une image virtuelle 
eiementaire de resolution 64. On divise alors chaque pixel en quatre de 
fa9on k obtenir quatre fois plus de pixels et k agrandir cette image virtuelle 
elementaire de fagon k la porter k la resolution 128. 

On projette ensuite dans le mSme plan les voxels du volume V V5 
(correspondsmit ^ la resolution 128), de fa9on k deflnir Timage eiementaire 
virtuelle de resolution 128 que Ton somme pixel k pixel avec Timage 
virtuelle agrandie precedente, egalement de resolution 128. 

On obtient alors une image sommee de resolution 128, dans 
laquelle on divise encore chaque pixel en quatre de fa9on k obtenir quatre 
fois plus de pixels, et de fa9Qn k obtenir done une image agrandie virtuelle 
de resolution 256. 

On projette les voxels du sous-ensemble VV7 de fa9on k obtenir 
rimage virtuelle eiementaire de resolution 256 que Ton sonune k Timage 
agrandie de resolution 256, de fa9on k obtenir Timage virtuelle 
intermediaire de resolution 256. 

On agrandit enfin cette image intermediaire de resolution 256, 
d'une fa9on analogue k ce qui vient d'fitre decrit, de fa9on k obtenir une 
image agrandie de resolution 512, que Ton somme pixel par pixel avec la 



2779853 



13 

projection des voxels du sous-ensemble VV3, de fa9on k obtenir en final 
Timage virtuelle IVi, de resolution 512, et correspondant k la projection 
des voxels du volume VV2 dans le plan de I'image acquise lAi. 

On effectue ensuite la difference entre Timage virtuelle IVi et 
5 Timage acquise lAi (dtape 48). 

On effectue enfin una r6tro-projection de cette difference (etape 
49) de faqon k corriger respectivement les valeurs de density des voxels 
des sous-ensembles respectifs VV3, VV5, VV6 et VV7. 

On obtient done des sous-ensembles corrig^s VV30, VV50, 
10 VV60etVV70. 

Ces operations de projections successives, de difference et de 
retro-projection sont effectuees pour chacune des images acquises lAi, 
On effectue ensuite une troisifeme iteration et Ton obtient en 
final une representation tridimensionnelle de Tobjet radiographie. 
15 L'algorithme de reconstruction d'image selon Tinvention opfere 

done en ce qui conceme les deuxi^me et troisi^me iterations sur environ 7 
millions de voxels qui se repartissent en environ 6 millions de voxels de 
resolution 512, 300.000 voxels de resolution 256, 82.000 voxels de 
resolution 128 et 100.000 voxels de resolution 64. 
20 On obtient par consequent une diminution du temps de 

reconstruction dans un facteur 4. 

Par ailleurs, le premier seuil de detection tj est predetermine en 
pratique de fa9on k avoir un nombre maximum de Tordre de 6 millions, 
pour les voxels de resolution 512. 
25 Ces voxels de resolution 512 correspondent aux vaisseaux 

sanguins radiographies, tandis que les autres voxels correspondent aux 
zones adjacentes et plus eioignees. 

En particulier, les voxels de resolution 64 correspondent, 
lorsqu'il s'agit d'une radiographie du cerveau, au demi-hemisph^re du 
30 cerveau non conceme par I'injection du produit de contraste. 
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REVENDICATIONS 

1 . Proc6d6 de reconstruction d'une image tridimensionnelle d'un 
objet k partir d'un ensemble d*images projet^es bidimensionnelles 
num6ris6es dudit objet obtenues par diff^rentes positions d'un appareil de 
prises de vues autour de Tobjet, comprenant un ^talonnage (30) de 
Tappareii dans lequei on elabore un volume virtuel entourant Tobjet et 
d6compos6 en voxels (Vi), une acquisition (31) dudit ensemble d'images 
projet6es bidimensionnelles num6ris6es, et une reconstruction (33) de 
rimage tridimensionnelle k partir des images bidimensionnelles projet^es 
acquises et d'un algorithme it^ratif de reconstruction alg^brique d'image, 
caract^rise par le fait qu'on effectue une premi&re iteration de 
Talgorithme avec une resolution d'image initiale pr€d6termin6e de fa9on k 
obtenir k Tissue de cette premiere iteration des premieres valeurs de 
densit^s pour les voxels dudit volume (VVl), on subdivise (44) une partie 
au moins (VV2) desdits voxels du volume virtuel en plusieurs ensembles 
(VV3, VV5, VV6, VV7) correspondant respectivement k des resolutions 
d'images diff^rentes multiples ou sous-multiples de la resolution initiale, 
et au sein de chaque iteration suivante de I'algorithme, on applique 
successivement ledit algorithme k partir de chacun desdits ensembles de 
voxels. 

2. Procede selon la revendication 1 , caract^rise par le fait qu'au 
cours de chaque iteration suivante de I'algorithme, on determine pour 
chaque image projetee acquise, une image virtuelle (IVi) en sommant (47) 
successivement des images virtuelles eiementaires respectivement 
obtenues par les projections dans le plan d'acquisition de ladite image 
projetee acquise des ensembles correspondants de voxels, 

3. Procede selon la revendication 2, caracterise par le fait qu'on 
eiabore les images virtuelles (IVi) eiementaires dans I'ordre croissant des 
valeurs de resolution en commen9ant par I'image virtuelle eiementaire 
correspondant k la valeur de resolution la plus f aible, et par le fait qu'aprfes 
avoir determine une image virtuelle eiementaire, on effectue une 
modification d'echelle de fa^on k obtenir une image virtuelle agrandie 
dont la resolution correspond k celle de la prochaine image virtuelle 
eiementaire k determiner, et on somme ladite image virtuelle agrandie 
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avec ladite prochaine image virtuelle 61^mentaire. 

4. Proc6d6 selon Tune des revendications pr6c6dentes, 
caract^ris^ par le fait qu'on acquiert des images projet6es 
bidimensionnelles ayant une resolution de base pr6d6temiin6e (r = 512), 
5 on choisit une resolution initiale (r/2) 6gale k un sous-multiple de la 
resolution de base, et on choisit parmi lesdits resolutions diff6rentes 
d'images, ladite resolution de base (r), la resolution initiale (r/2), et au 
moins une premifere resolution suppiementairc (r/4 = 128) sous-multiple 
de la resolution initiale (r/2 = 256). 

10 5. Proc6de selon Tune des revendications precedentes, 

caracterise par le fait qu'on eiabore un premier seuil de densite (t|) en 
fonction d*un critfere de selection predetermine, par le fait qu'on subdivise 
chaque voxel ayant une densite superieure ou egale au premier seuil en un 
premier nombre de voxels subdivises, le premier nombre etant defini 

15 compte tenu du rapport entre la resolution de base et la resolution initiale, 
tous les voxels subdivises formant un premier ensemble de voxels (VV3) 
correspondant k ladite resolution de base, par le fait qu'on regroupe 
certains au moins des voxels dont la densite est inferieure audit premier 
seuil et qui satisfont k un crit^re de regroupement predetermine, de fa9on k 

20 ^fonner des vox^^ ensemble un deuxifeme ensemble de 

voxels (VV5) correspondant k ladite premiere resolution suppiementaire, 
le nprnbre de voxels regroupes dans chaque groupe tenant compte du 
rapport entre la resolution initiale et ladite premifere resolution 
suppiementaire, et par le fait que les voxels dont la densite est inferieure 

25 audit premier seuil et qui ne satisfont pas au critfere de regroupement 
predetermine ferment un troisifeme ensemble de voxels (VV7) 
correspondant k ladite resolution initiale. 

6. Procede selon la revendication 5, caracterise par le fait qu'un 
voxel satisfait audit critere de regroupement si chacune des coordonnees 

30 du centre du voxel est un multiple de 2 et si la densite de chaque voxel 
voisin dudit voxel est inferieure audit premier seuil. 

7. Procede selon la revendication 5 ou 6, caracterise par le fait 
qu'on regroupe (46) les voxels regroupes qui satisfont audit critfere de 
regroupement de fa9on k former un quatri^me ensemble de voxels ( VV6) 

35 correspondant k une deuxi^me resolution suppiementaire sous-multiple 
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de la premifere resolution suppl6mentaire. 

8. Proc6d6 selon Tune des revendications pr6c6dentes, 
caract6ris6 par le fait que la partie (VV2) des voxels du volume virtuel, 
que Ton subdivise (44), est obtenue en 61iminant (41) les voxels situ6s 
dans une couche d*6paisseur pr€d6termin6e de la surface du volume. 

9. Proc6d6 selon Tune des revendications pr6c6dentes, 
caract^ris^ par le fait que la partie ( VV2) des voxels du volume virtuel, 
que Ton subdivise, est obtenue en 61iminant (42) les voxels isol^s dont la 
density est sup6rieure k un deuxi&me seuil pr6d6termin6, 

10. Proc6d6 selon Tune des revendications pr6c6dentes, 
caract6ris6 par le fait qu'^pr&s la premi&re it6ration de Talgorithme, on 
affecte k chaque voxel dont la valeur de density est sup6rieure k un 
troisifeme seuil pr6d6termin6 (t3), une valeur de density ^gale k la valeur 
maximale de density prise parmi les valeurs de density dudit voxel et des 
voxels voisins. 
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